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We analyze the prospect of measuring the neutron form factor of a nucleus through the detection 
of neutrino-nucleus coherent elastic scattering. We predict numbers of events in a liquid noble 
nuclear recoil detector at a stopped pion neutrino source. We discuss the precision required to 
distinguish between different theoretical models for the form factor. 
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I. INTRODUCTION 

Hadronic (pion, proton) scattering experiments have 
given us our first glimpse of the neutron distributions 
within nuclei. However, due to uncertainties and dis- 
agreements in theoretical fit calculations, the neutron 
radii of medium to heavy nuclei are still not sufficiently 
well knownp], Proton charge radii are accurately 
known from electron scattering [lj], but until we have a 
better understanding of neutron radii our picture of the 
nucleus will be incomplete. The neutron(proton) radius 
is related to the neutron(proton) density distribution in 
a nucleus. An accurate measurement of the neutron den- 
sity distribution is desirable not only for the sake of hav- 
ing a more complete picture of the nucleus, such knowl- 
edge will be applicable in other areas of physics as well. 
A precise measurement of the density of a nucleus will 
provide a better understanding of the saturation density 
of nuclear matter. Hence, a better knowledge of neutron 
densities could be useful for applications ranging from 
reducing uncertainties in atomic parity violation experi- 
ments to understanding properties of neutron stars (see 
[l], 0] and references therein.) 

In the expression for the cross section for a projectile 
scattering with a nucleus, the form factor gives a measure 
of the charge distribution in the nucleus. Specifically, the 
form factor is the Fourier transform of the density distri- 
bution of neutrons and protons. A measurement of the 
neutron form factor of the lead nucleus ( 208 Pb) using par- 
ity violating electron scattering has been discussedpj and 
an experiment tentatively scheduled at JLAB for 20080. 

In this paper, we present another way the neutron form 
factor of a nucleus could be measured — using neutrino- 
nucleus elastic scattering. The matrix element for the 
process of neutrino-nucleus elastic scattering^ is ob- 
tained by the superposition principle [5[: the individual 
amplitudes (with relative phase factors) for the neutrino 
to scatter off each nucleon are added and then the total 
is squared. For the case of a spin zero nucleus, ignor- 
ing radiative corrections, the differential cross section is 
given by Q 
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In this expression, Gf is the Fermi constant, E is the 
neutrino energy, T is the nuclear recoil energy, M is the 
mass of the nucleus, Qw — N — Z{\ — 4sin 2 6V) is the 
weak charge (sin 2 6W w 0.231), Q 2 = 2E 2 TM/ (E 2 - ET) 
is the squared momentum transfer and F(Q 2 ) is the form 
factor. Denoting the neutron and proton densities as 
Pn, P (r), the form factor is [7|] 
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From Eq. ((2]) we see that since (1 — dsin 2 ^) is small, 
a neutrino scattering elastically with a spin zero nucleus 
couples mostly to the neutron distribution. A measure- 
ment of the cross section for this process provides a mea- 
surement of the neutron form factor. Such a measure- 
ment would be complementary to the parity violating 
experiment with lead because it would provide additional 
data, obtained at different energy ranges and with differ- 
ent nuclei, that could be used to calibrate nuclear struc- 
ture calculations. 

To illustrate our idea, we consider the scenario of a nu- 
clear recoil detector at a stopped pion neutrino source. 
However, the idea is general and an experiment could be 
preformed at a low energy beta beamj^, as well. (Us- 
ing a nuclear recoil detector at a low energy beta beam 
was considered in References our exam- 

ple scenario, we choose the CLEAN detector[l3| and the 
decay at rest neutrino source at the Oak Ridge SNS [141 . 
Such an experimental setup was proposed in Ref. [l5j . 
Note that the idea of using neutrino-nucleus elastic scat- 
tering to measure the neutron form factor could have 
been conceived years ago, when the neutrino's existence 
and the neutral current were verified. However, the idea 
is practical in this age because neutrino beams and nu- 
clear recoil detectors are becoming a reality. 

The paper is organized as follows. In Sec. II, we discuss 
the neutrino flux at the Oak Ridge facility as well as the 
expected capabilities of the CLEAN detector. In Sec. Ill, 
we present the general formalism for calculating number 
of events in a detector and present events in the detector 
predicted by different analytic models for the form factor. 
In Sec. IV we analyze the potential for a detector to 
distinguish between form factors from different nuclear 
structure calculations. We give conclusions in Sec. V. 
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II. NEUTRINO SOURCE AND DETECTOR 

The spallation neutron source (SNS) at Oak Ridge Na- 
tional Lab will produce neutrons by firing a pulsed proton 
beam at a liquid mercury target. Among the fragments 
emitted will be pions. Negative pions will capture in the 
target, but positive pions will come to rest and decay 
into positive muons and muon neutrinos. The muons 
will also come to rest and decay, emitting a positron, 
electron neutrino and muon anti neutrino. The pion de- 
cay time is 26ns, and the muon decay time 2.2/^s. The 
neutrinos will be emitted isotropically from the target. 
The energy spectra of the v e and 9„ are precisely known 
since these neutrinos come from decay at rest; the 
will be mono-energetic. The nu-SNS (neutrinos at the 
SNS) proposal [lij calls for putting a neutrino detector 
in a room about 20 meters from the mercury target. The 
direction of the proton beam to the target defines the 
forward direction; the neutrino detector would be located 
to the rear of the target, perpendicular to the beam line. 
In the rear location, the background of neutrinos com- 
ing from any pion or muon decay in flight in the for- 
ward direction would be negligible. The pulsed nature of 
the beam will be advantageous in reducing other types 
of background. The time difference for the emission of 
v e 's and P M 's may also be used for distinguishing between 
events coming from these neutrinos and v^s as well as 
further reduction of background. An expected flux of 10 7 
neutrinos per sec per cm 2 of each flavor is expected at a 
distance of 20 meters from the target, where the neutrino 
detector will be placed. 

The CLEAN (cryogenic low energy astrophysics with 
noble gases) detector [lj] was conceived for the purpose 
of detecting dark matter particles and low energy neutri- 
nos. The detector will contain a liquid noble gas. When 
an elastic scattering event occurs in a liquid noble gas, 
ultraviolet scintillation light is given off. The CLEAN de- 
tector will use wavelength shifting film to convert ultra- 
violet light to visible light that can be detected by photo- 
multipiers. The detector will be able to detect neutrino- 
electron and neutrino- nucleus clastic scattering events. 
Some neutrino applications include measuring the p-p 
solar neutrino flux, measuring neutrino cross sections to 
search for neutrino magnetic moment contributions, and 
detecting neutrinos from a supernova. The detector con- 
cept is currently undergoing testing. A full sized detector 
is expected to have a mass of 10-100 tons of liquid neon 
or argon. 21 

An experiment to use a CLEAN detector at the SNS 
has recently been proposed |15j. The neutrino- nucleus 
elastic scattering process has never been observed, and 
so a first goal of this experiment would be to detect this 
process. Detection of this process could also be used 
to search for beyond standard model physics by mea- 
suring sin 2 #w, discovering or constraining non-standard 
neutrino interactions by looking for their contributions 
to the cross section, and looking for a neutrino magnetic 
moment contribution to the cross section. The potential 



of these physics applications was studied in Ref. [l5| and 
it was found that an experiment would be most sensitive 
to search/constrain non-standard neutrino interactions. 



III. FORMALISM TO CALCULATE EVENTS 

We now study the potential of an experiment using a 
CLEAN detector at the SNS to constrain nuclear physics 
theories. We demonstrate this by presenting calculated 
numbers of neutrino-nucleus elastic scattering events in a 
detector. Different nuclear models predict different neu- 
tron distributions, and hence different form factors for 
the cross section. Therefore, given the neutrino flux dis- 
tribution, different nuclear models will predict different 
numbers of events to occur in a detector. The questions 
we address are: how much do the event rate curves differ 
according to different nuclear structure calculation pre- 
dictions for the form factor, and what sensitivities must 
the experiment achieve to distinguish between these the- 
ories? 

We consider events of v e and Du scatterings. These 
neutrinos come from the muon decay at rest. Their spec- 
tra can be calculated from the expression for the differen- 
tial decay rate of the muon (see Ref. :17j) by integrating 
over the energies of the electron and v a for the v e spec- 
tra, and integrating over the energies of the electron and 
v e for the spectra. The normalized spectra are 
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where is the muon mass. These expressions give the 
probability that a neutrino is emitted with energy in the 
range (E, E + dE). Plots showing the spectra shape can 
be found in Ref. [15| . The number of events are obtained 
by folding the neutrino flux with the differential cross 
section. In particular, the number of events per nuclear 
recoil energy is 



m » /2 da 

f(E)—(E,T)dE, (5) 



Emin(T) 



dT 



where E min (T) = \{T + VW+WM) is the minimum 
energy a neutrino must have in order to be able to give 
the nucleus a recoil energy T, N t is the number of tar- 
get nuclei in the detector, and C is the total number of 
neutrinos per second per cm squared of a given flavor 
reaching the target. Multiplying the spectra f(E) by C 
gives the neutrino flux at the detector. 

We do not use form factors from nuclear structure cal- 
culations. Rather, we represent different density distri- 
butions using analytic expressions which come from Ref. 
[lH . The density distribution from Ref. [l8| models the 
nucleons as having a constant interior density and sur- 
face thickness s. The Fourier transform of the density 
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distribution gives the form factor [18j 

F(Q 2 ) = ^^e X p[-i(Q S f] 



(6) 



where Rq = R 2 — 5s 2 , R is the radius of the nucleus, and 
s is the surface thickness of the nucleus. The momentum 
transfer Q 2 has been given above. We are treating the 
proton and neutron distributions separately, as in Eq. 
@. The form factor we use is 



F(Q 2 
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Q IL 



[NF n (Q 2 )- Z(l-4sin 2 w )F p (Q 2 )} (7) 



were F n (Q 2 ) and F p (Q 2 ) are the neutron and proton from 
factors, respectively. We use the functional form of Eq. 
([6]) evaluated at the neutron RMS radius R n and the 
proton RMS radius R p for the neutron and proton form 
factors, respectively. 

We consider the case of a detector filled with argon iso- 
tope 40 Ar, with Z = 18 and N = 22. An experimentally 
determined value for the mean square charge radius of 
this isotope is < R 2 >~ 11.75 [19(. Therefore, we take 
the RMS value of the proton radius to be R p — 3.43. 
For each form factor we use the same value of s = 0.5 fm 
for the surface thickness of the density distribution. Of 
course, if the experiment is to be performed, one would 
use a specific nuclear structure calculation to self consis- 
tently model the nucleus and all parameters required to 
fit the data. The purpose of our calculation is to demon- 
strate the potential of doing such an experiment. It is 
sufficient to describe the nucleus with an analytic model 
and assign values for the parameters as we have done. 
(We have in fact modified the parameters slightly and 
found our results and conclusions do not change.) 

We are interested in using a measurement of the neu- 
tron form factor to distinguish between different nuclear 
structure calculations that predict different neutron dis- 
tributions. We are modeling the neutron distribution 
with an analytic expression. The parameter in this ex- 
pression that represents different theories is the neutron 
RMS radius. Thus we will predict what events are ex- 
pected in the detector, depending on what nature actu- 
ally has chosen for the neutron distribution in the 40 Ar 
isotope, by varying the neutron radius. 

A current figure of merit for the neutron radius is that 
it is known to within an uncertainty of 10% for medium 
to heavy nuclei, though there is debate in the literature 
on this issue, l|. For a nucleus with equal numbers of 
protons and neutrons, one may expect that the neutron 
and proton density distributions are equal and that R n « 
R p . However, because of Coulomb repulsion, the proton 
radius may be larger. For a large nucleus, where there are 
more neutrons than protons, we expect that the neutron 
radius will be larger. For our calculations, we first set 
R n = R p and then vary R p by ±20%, ±15% and ±10%. 
Even though neutron radii are likely known to 10%, it 
will be useful exercise for us to consider the cases of R n = 



R p ± 20% and i?„ 
section. 



Rp ± 15%, as will be seen in the next 



In Fig. (p} we show an example of events in the detec- 
tor. This plot shows the predicted total number of events 
that will occur from v e and D u elastic scattering, per keV, 
per year, per tonne of detector material, plotted against 
nuclear recoil energy in keV. (Note that tonne = metric 
ton.) The numbers for this plot were calculated for the 
case R n = R p . Events are plotted per tonne so that we 
may be general and not impose a size for the detector. 

Next we consider events for R n = R p ± 15% and com- 
pare these to events in the case R n = R p . In Fig. ©, 
events have been summed and placed in 10 keV bins; the 
first five bins are shown. The solid line shows total events 
in each bin for the case R n = R p , and the dashed lines 
show total events in the cases R n = R p ± 15%. The dif- 
ference in binned event rates for these cases is visible in 
the plot. We can see that for low momentum transfers, 
i.e. low recoil energy, the events converge. This is consis- 
tent with the normalization condition on the form factor, 
F(Q 2 = 0) = 1. (The matrix element of the vector cur- 
rent for a nuclear state, in the zero momentum transfer 
limit, is equal to the conserved charge associated with 
that current [13]. Since we have factored out the charges, 
see for example, Eq. J7|) and Eq. (fT]), the proper condi- 
tion for our expression in this limit is normalization to 
1.) The low energy threshold of the detector is expected 
to be about 10 keV, i.e. the detector is not expected to 
be able to see nuclear recoil events with energies below 
10 keV. When analyzing the potential of this experiment 
to distinguish nuclear structure calculations for the from 
factor, we will consider events above the 10 keV thresh- 
old. However, later we will comment on opportunities for 
distinguishing theories if the threshold of the detector can 
be made lower. 

As there are thousands of events per bin in Fig. © 
the scale makes it difficult to see how the events differ 
over the whole range of nuclear recoil energy. Therefore, 
we present Table U which shows for each bin range the 
number of predicted events for each of the cases R n = R p , 
R n = R p ± 20%, and R n = Rp ± 10%. The table also 
shows the percentage difference in events for the modified 
cases of R n compared to the case R n = R p . 



IV. ANALYSIS 

Using the table we can best gauge the potential for the 
experiment to be used to fit nuclear structure calculations 
for the form factor. The fit will depend on what the fi- 
nal error bars are for the data. Guided by the discussion 
in Ref. [HI], we consider the scenario of 10% system- 
atic uncertainty on the data. This can take into account 
uncertainties in the proton beam, neutrino flux reaching 
the detector, and the detector. Note that this systematic 
uncertainty need not apply to all energy bins. Events 
in some bins may be resolved better or worse. In our 
analysis we will consider some effects of backgrounds and 
detector efficiency. We remark here that backgrounds for 
some bins may be more or less severe, and the detector 
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might have better or worse efficiency for certain ener- 
gies. It should be kept in mind that the true potential 
of measuring the neutron form factor with this type of 
experiment can be evaluated once backgrounds, detec- 
tor efficiencies, and systematic uncertainties are better 
known. 

We now consider the case of R n = R p ±20%. From the 
table, we see that events in the bins with recoil energy 

> 70 kcV differ from events for the case R n = R p by more 
than 10%. Therefore, if the systematic uncertainty on the 
data in these bins is 10%, and there are enough events for 
low statistical uncertainty, then there would be potential 
for the data to be fit to nuclear models that predict the 
neutron radius to have values R n = R p ± 20%. Even 
if background cuts and detector efficiency would remove 
some of the events from being used, it is likely that there 
is still a large enough difference between events from the 
models that they could be distinguished. For example, 
in the 80 — 90 keV bin, there is a difference of about 50 
events between the cases R n = R p ± 20% and R n = R p . 
This should leave enough room for some adjustments for 
background cuts and efficiency. 

The exercise of considering the case R n = R p ± 20% 
is useful because it demonstrates that this type of exper- 
iment could be used to measure a neutron distribution. 
However, we considered a specific experimental setup for 
our example, and the capabilities and precision of this 
setup are still being determined. Obviously, if the ex- 
periment is more sensitive, a better measurement of the 
neutron distribution can be made. Therefore, we next 
consider what precision would be required to determine 
R n to better than 10%, the current figure of merit for 
how well the neutron radius is known. 

Consulting the table, we see that for the case R n = 
R p ± 10%, the events in the bins with recoil energies 

> 100 keV differ from the events for the case R n = R p by 
more than 10%. Therefore, if the systematic uncertainty 
on the data in these bins is 10%, then one might think 
the data can be used to distinguish nuclear structure cal- 
culations that predict the neutron radius to have values 
R n = R P ± 10%, from theories that predict R n « R p . 
However, in these energy bins, the models R n — i? p ±10% 
differ from the model R n ~ R p by only a few events. 
Since the detector will not be 100% efficient, and back- 
ground events will need to be accounted for, our calcula- 
tion indicates these bins will not provide useful informa- 
tion for fitting to theory. We therefore suggest a goal for 
the systematic uncertainty needed to measure the neu- 
tron form factor of 40 Ar with enough precision to extract 
R n to an accuracy better than 10%. With our example 
for an experimental setup and assuming enough events 
to have a low statistical uncertainty, the systematic un- 
certainty would have to be reduced to better than 5%. If 
this were the case, then data from bins > 50 keV could be 
used. For example, our calculation indicates a difference 
of about ±50 events between the cases R n — R p ± 10% 
and R n = R p in the bin 60 — 70 keV, with percentage dif- 
ference of > 6%. This suggests that if systematic uncer- 



tainties can be reduced to 5% or better, then even when 
background and detector efficiency is accounted for, the 
data can still be used extract Rn to better than 10%. 

Note that the calculated events we presented are given 
per year, per tonne of detector material. If data is taken 
longer, or the detector mass is larger than 1 tonne, more 
events will be available and hence statistical uncertainies 
will be reduced. Note also that other nuclei, for example 
Ne or Xe, could be used in the detector. The prospects 
of doing a neutron form factor measurement with other 
nuclei should be investigated. We remark again that our 
analysis was done for a particular example of an experi- 
mental setup. A different experimental configuration will 
obviously have a different potential for doing this mea- 
surement. 

In our analysis we have ignored contributions to the 
cross section from beyond standard model physics. How- 
ever, it is possible that new physics exists in nature. If 
this is the case then how could we know if a difference 
in measured events from predicted events stems from a 
modification to constants in the cross section or a mod- 
ification to the nuclear form factor? One feature at our 
disposal is that different form factors will have slightly 
different shapes, while a correction to say, sin 2 ^, would 
shift the cross section. A shift in the cross section would 
shift the number of events at all recoil energies, includ- 
ing the lowest recoil energy. However, a correction to the 
form factor would not effect the low recoil energy events. 
(Recall the form factor is normalized to one at zero mo- 
mentum transfer.) Such differences in the event rates 
could be observable if the threshold of the detector can 
be made lower. 



V. CONCLUSIONS 

We have considered a new idea for a method to mea- 
sure the neutron form factor of a nucleus — by detecting 
neutrino-nucleus clastic scattering. We illustrated this 
idea using a particular example of an experimental con- 
figuration — the neutrino source at the Oak Ridge SNS 
and a CLEAN detector filled with one tonne of liquid 
argon. We found that for this setup, it is possible to 
measure the neutron radius of the argon isotope 40 Ar 
to better than 10% if systematic uncertainties could be 
reduced to better than 5%. The idea is general in the 
sense that other nuclei could be used in the detector, or 
a different low energy neutrino source could be used. 
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Bin Range 
(keV) 


Rn — Rp 


Rn — Rp 

+20% 


% diff. 


Rn — Rp 

-20% 


% diff. 


Rn — Rp 

+10% 


% diff. 


Rn — Rp 

-10% 


% diff. 


0-10 


4756 


4714 


-1 


4793 


1 


4729 


-1 


4781 


1 


10-20 


3891 


3800 


-2 


3971 


2 


3832 


-2 


3946 


1 


20-30 


2884 


2772 


-4 


2984 


3 


2811 


-3 


2952 


2 


30-40 


2126 


2010 


-5 


2230 


5 


2050 


-4 


2196 


3 


40-50 


1549 


1441 


-7 


1647 


6 


1478 


-5 


1616 


4 


50-60 


1110 


1015 


-9 


1197 


8 


1048 


-6 


1169 


5 


60-70 


778 


700 


-10 


851 


9 


726 


-7 


827 


6 


70-80 


529 


468 


-12 


587 


11 


489 


-8 


568 


7 


80-90 


347 


301 


-13 


390 


13 


316 


-9 


376 


8 


90-100 


215 


184 


-15 


246 


14 


194 


-10 


236 


10 


100-110 


124 


104 


-16 


144 


16 


111 


-11 


137 


11 


110-120 


64 


52 


-18 


75 


18 


56 


-12 


71 


12 


120-130 


27 


22 


-19 


32 


19 


23 


-13 


30 


13 


130-140 


8 


6 


-21 


9 


21 


6 


-14 


9 


14 



TABLE I: Table shows binned events. All events are given in units of per (keV tonne year). Bins ranges are given in the first 
column in keV. The second column shows the number of events in each bin predicted by our model for the case R n = Rp- The 
third column shows the predicted number of events for the case R n = Rp + 20%. The fourth column shows the percentage 
difference between events for the case Rn — Rp + 20% and the case R„ = Rp- The fifth column shows the predicted number of 
events for the case Rn = -Rp — 20%; the sixth column shows percentage difference between events for this case and R n = Rp- The 
eighth and tenth column show events for the cases R n = -Rp ± 10%, while the ninth and eleventh columns show the percentage 
difference for the events of these cases and the case R„ — R p . All numbers are rounded to the nearest whole number. 



v e +v^ Scattering Events at SNS 




Nucleus Recoil Energy (keV) 



FIG. 1: Number of v e + v u events in detector in units of 
number per (keV year tonne) assuming 1 x 10 7 neutrinos per 
second of each flavor are emitted from the source. Neutron 
radius is set equal to proton radius. 
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FIG. 2: Number of v e + events per (keV year tonne) 
summed in 10 keV bins. Note the range: only the first five 
bins are shown and the verticle axis starts at 1000 events. 
The solid line shows events predicted for the case Rji — Ftp 
and the dashed lines show predicted events for the cases 
Rn =Rp± 15% 
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